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Abstract 



The antisymmetrization of the composite particles in cluster model calcula- 
tions manifests itself in Pauli forbidden states (ghost states), if one restricts 
oneself on undeformed cluster in the low energy region. The resonating group 
method and the generating coordinate method rely on a property of the norm 
kernel, which introduces some of the ghost states. The norm kernel has been 
usually calculated under the assumption that the inner wavefunctions have 
a simple Gaussian form. It is the first time that this assumption is tested 
by the rigorous way. In the 4 He+N system, we demonstrate a ghost state, 
which is calculated from a rigorous solution of Yakubovsky equations for the 
a particle. The ghost states calculated by rigorous and approximate methods 
turn out to have a very similar form. It is analytically proved that the trace 
of the norm kernel does not depend on the inner wavefunction we choose. 
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Since 1937 H, when the resonating group method (RGM) was established, it has been 
successfully applied to many light nuclei systems. The method is essentially based on the 
variational principle under the conditions that the clusters remain in the ground state in the 
low energy region, and the total wave function is totally antisymmetrized because of the Pauli 
exclusion principle. A typical example is the two a model of 8 Be. In 1970's the RGM had 
great successes {§] and the method was extended to the generating coordinate method (GCM) 
||, the orthogonality condition model (OCM) [[]], the fish-bone optical model (FBOM) J7|, 
etc. The Pauli exclusion principle plays an important role in the relative motion part of the 
wave function because it rules out part of the model space by an orthogonality condition. 
The Pauli forbidden states (ghost states) are generated by diagonalization of an integral 
kernel in the RGM. The integral kernel is known as the norm kernel (NK) M which is 
defined, for example, in the two cluster model as 

Af(r, ?) =< 0!0 2 5(f) |(1- A)\(f>KhS(?) > (1) 

where 0i, 02 and A are two inner cluster wavefunctions of the system and an antisymmetrizer 
for all nucleons, respectively, r is the relative motion coordinate. Namely, the ghost states 
u n (r) are eigenstates of Af with the eigenvalue j n = 1. 



dr'Af(r,r')u n (r') = j n u n (r) (2) 
The total wavefunction ^/ of the system 

\^>=A\4>i4>2X> (3) 

is orthogonal to the ghost states |0i02«n > 

< *|0i0 2 w n >=< <j>i<j)2X\A\(j)i(j)2U n >=< 0i02Xl(l - 1)|0i0 2 m„ >= 0, (4) 

where x is the relative wavefunction between the two clusters. 

Beyond 8 Be the a cluster model has been studied in 12 C f|-10], 16 [11|, etc. The correct 



treatment of the Pauli forbidden states is essential even in the case of bound states of clusters, 
where the neglect of the Pauli principle leads to an extreme overbinding. However, even with 
the condition they are over-bound, which is still a pending problem fr2HT4| in the a cluster 
model. 



It is analytically proved |Tj| that if the inner wavefunctions are simple products of Gaus- 
sian functions, then the eigenvectors u n of (0) become familiar harmonic oscillator functions. 
For the sake of simplicity four spinless cluster (a particle) system in 16 O has been studied 
i~T| using the Yakubovsky equations |16|]. Nowadays, it is possible to obtain rigorous solu- 



tions of the a particle wavefunction using realistic potentials ||17|| . Recent Progress of 4N 



scattering state is reviewed in jL8|. Therefore, it becomes possible to compare the ghost 
states which one obtains from the rigorous Yakubovsky solution of the a particle. 

In this letter we would like to choose the most simple case, the a + N system. 

The first 4 nucleons build up the ground state of the a particle, while the fifth particle 
is the spectator. The nucleons are identical particles, furthermore, the wavefunction <p a of 
a particle is normalized to < Q |0 a >= 1 and NK is defined as Eq. (P, 

M = 4 < Q n 5|P 45 |0a0„5 > • (5) 
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P45 is the particle exchange operator (4 and 5). Fig.|I] suggests the picture of three cluster 
system (3N + N +N). The calculation of the NK is done with Jacobi coordinates which we 
show in Fig. 0. The relative Jacobi momenta are prepared and the relations between the 
Jacobi coordinates are 



Pi 



<12 + jqi, 



P2 



^ 1 ^ 

Qi + 4^2 



(6) 



where Pi and (i = 1, 2) are Jacobi momenta of 3N+N relative motion and 4N + N one, 
respectively. 

Our way of calculating the NK is very similar to the calculation of a leading (Born) term 

For example, in the text [pUl 



of the Alt-Grassberger-Sandhas equations |T9|J20[| . For example, in the text [^JJ they treat 
the Born term Z by the partial wave representation with a function F c . In our calculation 
it is simply replaced as 





/■OO P 


j d cos 6 


I dx j 




Jo Jo 



dy<f) a:Nl {x,y, \pi\)(j> a -.N 2 {x,y, |p 2 |) 



P c (cos9) (7) 



where the angle 9 is between the vectors <fi and g*2, N\ and N2 the state channels of the 
partial wave, and Pc the Legendre function, x and y are rests of Jacobi momenta which 
describe the motion of particles (1,2 and 3) inside of the a particle. The numerator 4 derives 
from (|5|). 

As an example, the a wavefunction |l7j of the Argonne potential (AV14) [^2| is applied, 
and we take the case of total spin J = l/2 + . For the sake of simplicity we assume the spin j 
of 3N is almost l/2 + (in fact, 94.9% for the case of AV14 potential ), therefore, the angular 
momentum between clusters a and neutron is S-wave. This leads to £=0. 

Under this choice of the partial waves the recoupling coefficient N2 
one gets 



21 is 1/4 and 



■A/o (91,92) 



dcos6 



dx / dy(j) a . [1/2 +](x,y,pi)(p a . [1/ 2+](x,y,p2) 



d cos 9M(pi,p 2 ) (8) 



where the subscript "0" of the norm kernel means the angular momentum of q and the kernel 
Af will be used later. 

The ghost state is shown in Fig.|[ For our calculation the eigen value 70 of Eq. (0) is 
not exactly equal to one, but 0.937 (if it is renormalized by the abovementioned 94.9%, 70= 
0.987). The solid line is the ghost state calculated from our Yakubovsky solution <p a , 
comparing to the dashed line from usual Gaussian function, 



128 ul N 



1/4 



exp(-u; Q ivr 



(9) 



with uj a N = Q x (4 • l)/(4 + 1) where Q is a common shell model mode (0.275 fm ~ 2 ) ||. 
They are normalized by J °° u\(r)r 2 dr = 1. In the short range our ghost state is smaller 
than Uq. The repulsive core of realistic potentials reflects in this range. This behavior is 
similar to that of correlation functions [p3 |. Beyond 4 fm our is bigger than Uq because 
in general the Gaussian function is more quickly decreasing than exponential one. We also 
show them in the momentum space (see Fig||). Here the repulsive core manifests itself by 
a node at w 2 fm ~ l which is absent in Uq. Overall they agree well. 
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To find the most realistic width parameter Q we optimize R = \ < u^\uq{Q} > | 2 x 100 
[%} in Fig.(||). We could recommend Q = 0.24[fm~ 2 ] of the Gaussian width parameter which 
is similar to ft=0.275 [fm~ 2 ] §. 

In table | we summerize the biggest eigenvalues in Eq. (H). Analytically we find in the 
Gaussian case eigenvalues 7 n = (— 4)~( 2n +°) , n = 0, 1, 2, . . . The realistic NK has got 
a similar spectrum. We compare the states u\ and uf for n=l in Fig. || It is remarkable 
that in this case the realistic ghost state has more structure though the eigenvalues are very 
similar. 

The matrix traces Tr [Ao] are given 



G 



Tr\j\ Q 
Tr[Afo] Y 



f 1 

Af (q,q)q 2 dq = 1.0125 



n=0 

roo 



16 
15 



1.0666 



(10) 



If the wavefunction of a particle is renormalized by only j = l/2 + , / °° p 2 dpAf(p,p) = 1 ( 
0.949 : original norm ) we get Tr[A/o] y = 1.0666 which must exactly be the number of the 
Gaussian form. Because it is analytically proved that the trace Tr [Ao] does not depend 
what kinds of the a wave function we choose: 



Tr 



- / dcos9J\f(\ o, 

2 7-i V V16 2 y ' 



17 cos^ , 

Ye + — q) 



dcos8- 



17 I cos t 
16 2 



16 

15' 



q 2 dq 



(11) 



We illustrate both NKs (Fig. [7| for Af^ and Fig. || for A/" y ). The Gaussian case is 
analytically given, 



A" (<?,<?') = — 1 
qq' 



exp 



17 
48fi 



(? 2 + ?' 2 ) 



sinh 



(12) 



The shape is so similar that the difference (A/" y — A/" G ) is also shown in Fig. |9|. 

Although there is only a single ghost state in a-N system, in general, the cluster-cluster 
effective interaction in light nuclei has a lot of ghost states. In this simple case we could 
find some remarkable differences in the eigen state (n=l) and the eigenvalue for n=2 which 
might effect RGM calculations of systems with A > 5. For a most probable case such a 
Pauli blocking will be applied to the a-n-n three-body model system. There are already 
some applications by using some Pauli methods. 

It will be important benchmark calculations for more nucleons system to look into the 



ghost states using rigorous solutions [26] from Few-Body Physics. 

Note that here we discuss the Pauli forbidden state which is different from the spurious 
state of the Faddeev calculations 



The naming of spurious state has been used a lot 
in many places, even if a cluster model has no inner structure the ghost states appear in the 
model and they are interpreted kinds of spurious states. We should not confuse spurious 
states caused from Faddeev decomposition |27|-p9l. In this paper we simply take the physical 
Yakubovsky solution of a particle to test quantitatively how precise the former Gaussian 
norm kernel is. 
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FIGURES 



FIG. 1. The diagram of the norm kernel. 
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FIG. 4. The ghost state (n = 0) in momentum space. The solid (dashed) line is Uq (q) (uq 
The disconnection of the solid line causes from change of sign. 



100 




1.22 n 0.^ 
D, [ fm- 2 ] 

FIG. 5. The percentage R of the ghost state as a function of oscillator parameter Q. 
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FIG. 8. Norm kernel A/Jf. 




FIG. 9. The difference of the norm kernels(A/o^ — 
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TABLES 



TABLE I. Eigen values the norm kernel. 



n 


In 


iln 1 ) oiu n 


In 


( 7 ^)of«G H 





0.937 


(1.068) 


1.00000 


(1) 


1 


0.0663 


(15.09) 


0.06666 


(16) 


2 


0.00753 


(132.) 


0.00391 


(256) 
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